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ABSTRACT 


The direct current resonance probe technique was applied to a 
moderately high density positive column xenon discharge in the presence 
of moving striations. The resonance frequencies obtained were of the 
order of 500 megacycles per second, corresponding to electron plasma 
frequencies of the order of one to one and a half gigacycles per second. 
The results of the experiment indicate that the resonance probe tech- 
nique may be a useful method for the determination of electron tempera - 
ture and density, but not electron-neutral collision frequency, in 
moderately high density plasmas, Experimental difficulties are dis- 
cussed, and recommendations for the course of future work with 


resonance probes dre presented: 
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1. Introduction. 

The resonance probe technique was first developed in 1960 by 
Takayama, Ikegami, and Miyazakif[l], as a means of quickly deter- 
mining electron temperature, density, and collision irequencey in a 
plasma. 

The method consists basically of applying a high-frequency signal 
to a Langmuir probe and noting the behavior of the collected direct 
current. Referring to Figure 1, the direct current shows threceance 
tinct regions as the frequency of the applied signal is varied: the first 
is an increase in the dc level in the region below the resonance fre- 
quency, independent of the signal frequency, from which the electron 
temperature can be determined; a noticeable peak at the resonance 
frequency; and at higher frequencies, a rapid fall-off of the direct 
current to the level corresponding to no applied signal. From the 
position of the resonance frequency, electron plasma frequency can be 
calculated, and from the half-width of the peak, the electron-neutral 
collision frequency can be determined directly. The basic circuit for 
the dc resonance probe method is illustrated in Figure 2, 

In this experiment, the resonance probe technique was applied to 
a moderately high density (on the order of 10° electrons as xenon 
positive column with pressures on the order of 10 COrr 

There were major differences in the type of experimental conditions 


here as compared to most previous work. Most previous work has been 
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done in '"'back diffusion'’ plasmas, with electron densities on the order 

6 -3 
of 10 cm _., andthe type of gas has usually been mercury or cesium 
vapor. In the present investigations, in addition to the high density 
xenon, there were moving striations present in all ranges covered. These 
striations had wavelengths of seven or eight cm and frequencies of 
25-100 kilocycles per second, with no apparent linear dependence on 


plasma parameters, They are ionization type waves which usually move 


from the anode to the cathode through the positive column [13]. 
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2. Resonance Probe Theory. 
Takayama and Ikegami [2] believed that the resonance frequency 
occurred at the electron plasma frequency, which is a function of electron 


density only: 


2 Lye 
(1) ae ( Aen = 22 8.97x10> (atesuai ny ee 


where a is the electron plasma frequency, n the electron densityaimot 
necessarily constant), and € and m the electron charpe and mass, 
respectively. Subsequent investigations [5] [7] [8] [9], however, have 
shown that the dc resonance peak always occurs below w 

For superposed signals of form A: sinwt, below the resonance 
frequency, the increase in direct current 1s explained as followec wane 
electron current density, ij , to the Langmuir probe, is given by the 


well-known relation: 


(2) ic 2 n( Ko Wye KT 


aie = 


where K is the Boltzmann constant, T the electron temperature, and 
V is the potential difference between the plasma and the probe. The 
subscript a will’reter to the current densivy duc tone prose bias, Vv, 
which is held constant in resonance probe investigations. 

The electron velocity distribution function is also assumed to be 


isOtropre, and is determined by; 


ait 





(3) f(y) -= Nr ii 


mere vis the electron velocity. 
According to Ikegami and Takayama [c], the direct Current increases 
67> 1s then given by the second derivative of the electron current densit, 


Versus pPOtLential driterenmce curve, 


(4) oe 


If the distribution is assumed Maxwellian, the dc increase can be 
Written: 

eA 
5 i as I sas - | : 
(5) SU | Gee ) 
where I, (2) is the zero order modified Bessel function, and A is the 
amplitude of the superposed signal. From this the electron temperature 


Canene Celermined,: 
(6) i (= -1 = 7 
O 

For frequency ranges above the resonance frequency, the decrease 
in 6j to zero is explained by Ichikawa [ 3] as being due to the effect of 
the oscillating potential cancelling out as an average, since the electrons 
cannot follow the fluctuating electric field, the transit time, Ue ~of thie 
electrons across the probe sheath, is much lomger than thegperiod of 


the potential oscillation, w 


Solving the equation of motion for the electrons: 


v4 


(7) mx =- eE + eE, sin(wt-a), 


A 
where E- = a : EB, oe is the sheath thickness and q@ is a phase 
O 
angle due to the plasma impedance defined by tana = € _ / ec. where 


.. and oF are the real and imaginary components of the plasma permit- 
tivity, respectively; and x is the distance of (he elect nonsmerom tie probe, 
Considering an electron initially at the edge of the sheath, which will 
require a minimum velocity a (hereafter called the "eritical velocity”) 
to reach the probe, leads to the following initial and boundary conditions 
for the solution of the equation of motion: 

<= seo = =Vrat t -.0ewand <= Oat x= 0, 


Men thevcritical velocityis derived: 








(8) v ee © ae (cos a - cos(wt, (a) - a) ). 
C mo o 
Then: 
(9) — saath fo v. (a, 0) vi(v) dvda. 
O O 


Solving equation (8) for the critical velocity in the low frequency 


region yields: 


Ze v2 


(10) v, (a) = (= (V- Asine ) ) 


Substitution of equation (10) into equation (8) leads to equation (5). 


Integrating equation (7), one obtains: 
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e E eE ek 
o 62 l l 
fs o- t cos (wt -a) - sin @ 
mw 











+ sin (wt - a) ) = Le 


limthevniphetrequency limit. the trancit time) Of clectrons across the 


sheath is found from equation (11): 


; 2ms 1/2 
(12) T= lim (t) = ¢ = 
ww) —+> © O 





Substituting equation (12) into equation (8), yields: 


i ore eA Zevaelpa 
Fo) CE ea ee a ae pet ey =~) 


eA l 
“ee kT 'Y) —). 





from which §6j approaches zero as ws becomes greater than v. 

ln thepireduency region at and near time mesonance frequency, acne 
mechanism is not as straightforward, and no universally accepted theory 
is available at present. A summary of some of the existing theories 
which are applicable follows. 

The plasma is assumed to be homogeneous, and for simplicity in 
development, plane parallel probe geometry is used. The plasma 
system can be represented schematically as in Figure 3(a), with the 
high frequency potential applied across the probe and a reference elect- 
rode, usually the cathode. The resulting equivalent circuit is illus- 
trated in Figure 3{b), after Levitski and Shashurin [4], Uramoto, et al 
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[5], and Mayer [6]. 

It has been confirmed that both the rf and dec current components 
show a maximum below the electron plasma frequency, and the the rf 
current drops to minimum near the electron plasma frequency. The 
former case (maximum) is interpreted as a series resonance, and the 
latter as a parallel resonance. Uramoto, et al [5], showed that the 
rf component increased monotonically at frequencies above the electron 
plasma frequency, while the dc component, 6j, as previously discussed, 
drops rapidly to zero. 

As the present investigations were concerned primarily with the 
dc resonance, the development will follow Ichikawa [3] and Ikegami 
and Takayama [2], realizing that the theory must be applied to the 
probe-plasma series capacitive circuit, and that the resonance fre- 
quency, w ae is not the same as the electron plasma frequency, 

If an external oscillating potential is applied at a frequency near 
the electron plasma frequency, the much slower ions present in the 
sheath about the biased probe will not be able to follow the oscillating 
electric field and their distribution will be a function only of the probe . 
bias, V. The electrons in the sheath, however, will follow the applied 
oscillating electric field and the resulting charge fluctuation in the 
sheath acts as an oscillating dielectric layer. The external perturbing 


field, oe , Calmoewexpressca as: 


(14) B act 6 ta ~ D(x) sin (wt-a ) , 
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where 


dk , 


(15) Dyed ee 4 [ eo sinks oikx 


k 


= © 


where k is the wave number. When these latter equations are used to 


solve the linearized Boltzmann equation: 


of af 
Ga a a On 
at ax m av m ext av at Eoll.., 


where f is the distribution function in an equilibrium state and f is 


the perturbation due to the applied signal, and also 


(17) o=- ane fia, 
Ox 


the following solution is obtained: 


af , 
(18) f(x, v,t) = — — af dw ABE Qe { a(w t+w)e. 
m av Oo Oo 
1a 
- a (WO -w)e i 
cael ge ikx ] 1 sin ks 
at w -kvtiv é(k, w_) ks 
O Cc O 


where WS is taken as a dummy variable, and a is the effective efectron 
collision frequency. 
When the mean path is much greater than the sheath thickness, the 


integral over k in the last equation can be estimated asymptotically: 


] ikx 1 l Sin ks l 1 
1 ie — dk Sn. ee ee oe eg ee 
7) on | ~  o-kvtiv. e(k,w) ks wtiv e(o,w ) 
s27o Cc O Cc Oo 


iG 


= 
=. 


Due to the time dependence of the perturbing rf signal, equation (9) is 
evaluated at t = 0, and assuming Maxwellian distribution of velocities, 


the total distribution function becomes: 


(vj) = f(y) i (x, vet = 0) 





The real and imaginary components of the dielectric permittivity are 


defined by: > 
" p 
(21 a) Pe) | ae 5 5 
yp =v 
Cc 
and 
Z 
Ww cwv. 
= p 
fel») aw, De 2 
W - V Ww - Vv 
Cc Cc 


Substituting equation (20) into equation (9) yields the final result: 


ey, eV \ > 
Ka Te, pea 
ss 2 


sone oe 
( ) J, + Sil) ip Me T(w) t+ io ae 





2 
oO +r Vv 
Cc 


a 
ot 


where: 
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+ € w (€ +6, ) 
r i o'r i 
°r Ve 7 
oth Fhe) ew) 
é 2 2  “ € ¢ =f 6° 
1 i r i 
where MD issune Depyveslenorh,. aid: 
RI Tee 
~ x 
= Aa Cc 
(23) 1) { e da 
O 
, ve eZ 
where the dimensionless critical velocity, x (a, Or) = (ey Ke) v (a, w), 
and 
s ( 1 20 a idee Si 
(24) ( wb) ) = | \~ xe. * dx( a » ae. 
c(w) 2 II cos @ 
O x (a, w ) 


Equation (22) is a general solution for total de collected by the 
probe, Using the critical velocity for the low frequency limit given by 
equation (10), the first term of equation (22) reduces to equation (5), 


If the approximation 


Zz Z 
2 = - 
(25) io Oe Ne sR ei at 
c 
x 
Cc 
is taken, equation (22) reduces to: ° 
Z 
IDE eee 5 Vea , De eA, eA 
JF dy 2 SW (s ecu 1: Ga 


Equation (26) can also be derived by solving the fundamental equations 
for linear response of the plasma electrons, equations (16) and (17), 


directly for dc current density, assuming that 


eV/KT is much greater than one, that A/V is much less than one, and 
Z 
that ( Ue /w) is also much less than one. 


Furthermore, if equation (26) is solved for the resonant peak 





herent, Tee assuming w = W = oe : 
W IN 
(27) — ae D é. eA ; (3) 
Jy oe YX Pea he im 


where IL, (z) is the first order modified Bessel function. 
The heli width ef the resonance peak determines the ellective 
electron-neutral collision frequency, the relation as shown by 


Wakayatmenmet al [-2il|, 1S: 
(28) Aw = @v 


This has been experimentally verified to within an order of magnitude 


by Cairns [7], using a low density mercury vapor plasma. 


ig 


3, Relationship of Resonance Frequency to Electron Plasma Frequency. 
As previously noted, the resonance peak of rectified current always 
occurs below the electron plasma frequency, and that this is due to 
series representation of the sheath-plasma interaction. 
Mayer [6] developed a theory for plane probe geometry and cal- 


culates the complex permittivity - of the plasma capacitor as: 


fhe 
2 2 
” 2 W ae fe 
€ - 6 € = = 
ey) eff a i Sy se. "'p ) Zz Z 22 
wW-w )+tvw 
r Cc 
where 
Z 2s Z 
a) * = pets “p 


This model has been experimentally confirmed by vonGierke, et al [ 8] 
for a low density cesium vapor plasma with unequal sheath thicknesses 
Alster properamd nererence electiode { S, and Sy: respectively) . 


Harp and Crawford [9] expanded this theory to spherical geometry 


and calculated: 


T 
= 


esd) Ww 





where R is the spherical probe radius and a Mn Heplaces Ss, and) cals 


a constant of proportionality. 
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4. Experimental Results and Conclusions. 

The resonance probe technique was applied to a xenon positive 
column discharge in a 60 cm Pyrex tube of 0.8 cm inside diameter. 

The probe was a 0.195 mm radius right-angle cylindrical tungsten probe 
located halfway between an oxide coated cathode and a tungsten rod 
anode. Neutral gas pressures of one to four microns of mercury were 
used in most of the experiments with electron densities in the plasma of 
the order of to” Beas . Tube currents were 20, 30, and 50 milli- 
amperes. 

Early attempts to find a consistent resonant increase in collected 
@irect current were unsuccessful.” it was believed that, in accordance 
with Harp and Crawford [9], the probe radius was too small in relation 
to Mp? thereby damping out any resonant peak. In order to eliminate 
this problem, a second probe was installed in the system near the right- 
angle probe, towards the anode. This was a thin tantalum band, two mm 
wide, which was mounted around the inner wall of the tube. The 
apparatus is diagrammed in Figure 4. After installation of this band, 
resonance peaks were observed, both on it, and on the original cylindrical 
probe, when either probe was at floating potential and the other biased 
near floating potential and with an applied rf signal. It would appear 
that the presence of the second conducting surface is essential to ob- 


servance of resonance frequencies in this type of discharge in that it 


provides a closer reference electrode than either the anode or cathode. 


fag | 


The radial profile of the electron density is assumed to be parabolic, 
and that the electron density at the wall is approximately 60% of that at 
taescencer of the tube, A typical plot of 6j/j,, versus the applied 
frequency is shown in Figure 5, with results obtained from it listed in 
aple\; 

It was expected that density variations caused by moving striations 
would broaden the resonance peak. This was confirmed in all experi- 
mental conditions as illustrated in Table 2. For calculation of v it the 
electron temperature as obtained from d(ln j,) / dV plots was inserted 


into 


ee Pe Nee 
a Ree Ge: iene! 


WheGespolssthe neutral gas pressure, in torr, and P . is the probability 
of collision at one torr, [10], uis the average electron velocity, and 
h E is the clectron mean free path: 

The expected peak height, bj. ie , as calculated from equation (27), 
did not show a consistent relation to those measured, but, instead, 
varied randomly. It was confirmed, however, that the magnitude of the 
peak height, and the de increase observed below the resonance frequency, 
were proportional to the amplitude of the applied signal. The results 
of peak height comparisons are given in Table 3, the amplitude of the 
applied signal was 0.7 volts. 


Thesresultssef tne measurement of the resonance irequency wien 


ape 


respect to the calculated electron plasma frequency appear to be closer 
to the model of Mayer [6] than to that of Crawford and Harp [9]. As 
shown by vonGierke, et al [8], the sheath thickness is assumed to be 
a function only of the Debye length, Mp? IRC 55 Siaeaay 
For the right-angle cylindrical probe (CP), the value of a is 2. 66 


for closest fit of the ratio of measured w tom@alculated Wi? using 


ae 


ve )/dV plots. For the wall, or 


equation (30) and Be obtained from d( 
mine. probe (RP)y the value ot ais’ 5.5. .in the latter case mips 

P 
approximated by determination of the random ion current density, dis : 


from Langmuir probe curves, and the relation from Glasstone and 


Bovbero (lilies: 


(33) n. ~ j.. (0.40 aja =< me ; 


1 A aE mM; 
1 


where m., is the ionmass. This assumes charge neutrality, i.e., the 
ion density, n =n. 

A sample of results of the above measurements is given in Table 4, 
and the ratio of calculated w to measured W is plotted against experi- 
mentally measured w in Figure 6. 

In the frequency region below the resonance frequency very poor 
agreement with theory was obtained. The increase in the dc level ap- 
peared to be independent of applied signal frequency, but a large random 
scatter was noticed, and, at best, only an estimate of 8/5, could be 


made. This was most probably due to drift of i and to the difficulty in 


Zo 


Beciratery 1eading the changes in the total current density, j, since 
6j=j- J was very small. Sample values are listed in Table 5, 

In the region above the resonance frequency, §j did not drop rapidly 
to zero, and appeared to have one or more resonance peaks, These 
peaks, as suggested by Ikezi and Takayama [12], may be Tonks - 
Dattner resonances, 

Since a resonant peak can be detected even at these comparatively 
high densities, it is concluded that the method of resonance probes may 
be useful as a diagnostic technique in positive column discharges for 
determination of electron temperatures and densities, However, due to 
the presence of moving striations, it is believed that the method may 


vyielduan average Collision irequcene von (ncuplacimd, 
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5. Recommendations for Future Investigations. 

In future work it 1s recomimended that ther! current be examuned 
with the ultimate goal of developing a simple method of direct determi- 
nation of electron plasma frequency. 

It is,further recommended that the de resonance pRabe acain be 
attempted in the absence of a second close electrode, “Ine problem ot 
electron temperature determination using the resonance paobe technique 
may be overcome if a leveled-power rf generator is obtained, anda 
method can be devised to eliminate the drift of Jo , and to measure only 
6j rather than j. The apparent resonant peaks above w should be more 
carefully investigated to determine whether or not they do actually 


represent Tonks - Dattner resonances or some other phenomenon, 


fags. 


Or 


ie 
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FIGURE 1 


Theoretical Plot of the Increase of dc Current 


Versus Frequency of the Applied rf Signal 
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FIGURE 2 


Basic Circuit for the Resonance Probe 
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(a) The Model Proposed for the Plane Probe 
Plasma-Sheath System 


(b) The Corresponding Equivalent Circuit 
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Analysis of Figure 5. Basic parameters are: xenon positive 


column discharge, pressure of 2.15 microns of mercury, and tube 


current 20 milliamperes. 


Resonance frequency, wo (ee i: 


Half width, Aw/2Tl: 


Collision frequency, v i i245: 


Resonance peak height, 65/4, 


Electron temperature, KT/e: 


490 megacycles per second 
100 megacycles per second 

50 megacycles per second 
0.090 


5.05+ 1. 50 ev 


From the corresponding j versus V probe curve, the following compara- 


tive data is obtained: 


Eke ctron density, n; 


Oy LOX io ae 


Electron plasma frequency, w /2: 1320 megacycles per second 
P 


Collision frequency: 
Plectron temperature - 
Blectron ebye length, ji D: 


Probe petential, V: 


32.1 megacycles per second 
4.05 ev 
OF OE segh on! 


23 volts 
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Representative Values of Measured Resonance 
Peak Height and Comparison with Calculated Values 


Ta a ae 
Electron Density Eee ae j,/dg) 


3 Ratio of 
(cm ~) Measured (1) Calculated (2) Column 2 to 1 
Cylindrical Probe 

10 

Cuoe x10 0.046 0. b22 2, 65 
10 

Zate x) 0.090 OLFnOS l. i 
10 

2.18x10 0.119 0.098 0. $2 
10 

3.70x10 0.078 0.158 Z. 03 
10 

4,.36x10 0.118 0.046 0. 39 

Ring Probe 

10 

1.03x10 0.038 0.0305 0. 80 
10 

ie 0 0.081 O40S 75 0. 46 
10 

2,01 < 10 0053 0.0732 1, 368 
10 

2,07 x10 0.081 0.0425 0. 525 
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PAB 5 
Representative Values of Electron Temperature 


as Determined by the Low Frequency Signal DC Increase 
Compared with Corresponding Values from Langmuir Probe Curves 


Electron Temperatures 


Pressure 
Current (microns of Resonance Langmuir 
(ma) mercury) Probe Probe 


Cylindrical Probe 


a0 Nes 9.82 4, 30 
20 Zale 0S 4.05 
20 2.4 35405 3. 80 
20 4.5 Seen Sao? 
Ring Probe 
20 le 2: (ase 4, 74 
20 S20 3, 54 4, 34 
20 a2 6. 43 4, 34 
50 Fa 8, S298 Bye le)S 
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